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Abstract. The cosmic baryon budget is revisited using modern obser- 
vations that have become available since our first publication. I also 
present an estimate for the heavy element abundance. An increased ac- 
curacy in the accounting of the baryon budget reveals 'missing baryons', 
which amount to ~ 35% of the total. This would provide an interesting 
test for models of the cosmic structure formation. 



1. Introduction 

The evolution of the Universe and the formation of cosmic structure redistribute 
dark matter and baryons. Thus the present day distribution of the matter 
reflects the history of the Universe, and gives us information as to how the cosmic 
structure formed. For this reason the matter distribution is taken as a useful 
constraint for models of formation of large-scale structure and galaxies. We have 
published an accounting of the cosmic baryon budget in 1998 (Fukugita, Hogan 
& Peebles 1998, hereafter FHP). Since then, much progress has been made in 
relevant observations, which include the Sloan Digital Sky Survey (SDSS: York 
et al. 2000), Wilkinson Microwave Anisotropy Probe (WMAP: Bennett et al. 
2003), several HI surveys, and others. In this talk, I attempt to update the 
present-day baryon budget using modern data, and discuss some issues relevant 
to cosmic structure formation. We write the Hubble constant as Hq = WOh km 
s _1 Mpc -3 , but adopt Hq = 72 km s~ 1 Mpc~ 1 when h is not explicitly denoted. 



2. Baryons in stars 

The basic data used to estimate baryons in stars are the luminosity function 
(LF) and the stellar mass to light ratio {M s /L) of galaxies. The most accurate 
LF was derived from the SDSS for five colour bands (Blanton et al. 2001; 
2003; Yasuda et al. 2003). The first LF from the SDSS given by Blanton et al. 
(2001) is based on earlier data for the northern equatorial stripe of approximately 
200 square degrees, giving the global luminosity density C r = (2.58 ± 0.28) x 
10 8 /iLq(Mpc) -3 . It turned out, however, that the surface density of galaxies in 
this region for r < 17.9 mag is somewhat overdense compared to the mean. The 
First Data Release of the SDSS (Abazajian et al. 2003) now covers 2200 square 
degrees, and the LF derived from these data gives a somewhat smaller value 
C T = (2.32 ± 0.25) x 10 8 /iL o (Mpc)~ 3 (Yasuda et al. 2003; see also Blanton et 
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al. 2003) in the r band. The corresponding luminosity density in the z band is 
C z = (3.9 ± 0.6) x 10 8 /iL o (Mpc)- 3 . 

The most extensive analysis for the stellar mass to light ratio is that by 
Kauffmann et al. (2003) using 10 5 SDSS galaxies. The estimate of the stellar 
mass using a population synthesis model depends on metallicity, age, and the star 
formation history, and also on the initial mass function (IMF). They estimated 
the probability distribution of each parameter in terms of the Bayesian analysis 
using five colour photometric data. In our consideration we limit to bright (r < 
15.9) galaxies, the median redshift of which is 0.05. The value of M/L z ~ 1.85 
for luminous galaxies with M z < M* —0.8, and it decreases gradually to 0.65 for 
galaxies with M z ~ M* + 3. The LF-weighted mean is (M/L z ) ~ 1.5 to which 
a 20% error is alloted. 

There is still a significant uncertainty in the IMF, especially for the subsolar 
mass. Kauffmann et al. assumed the IMF of Kroupa et al. (2001). Another 
typical IMF is that of Reid et al. (1999) from DENIS and 2MASS surveys. 
The former shows an increase of the IMF down to M = 0.1M©, whereas the 
latter flattens at M < IMq. This leads to a 40% difference in the integrated 
mass. We take the geometric mean of the two as our central value, allowing 
for ±20% errors. The Kennicutt (1983) IMF gives the integrated mass close 
to our adopted value. The IMF for < O.IMq is even more uncertain, but the 
contribution from this region is small (< 4%), provided that the IMF declines 
towards smaller masses (e.g., Burgasser et al. 2003). The value adopted here is 
1.10 times that used in FHP, which employed the subsolar mass IMF of Gould, 
Bahcall & Flynn (1996). The Salpeter IMF, when cut off at O.1M , gives the 
integrated mass 1.50 times the adopted value. 

Assembling these inputs we obtain 

O star = 0.0025 ± 0.0008, (1) 

which includes dead stars. This value is compared to 0.0019—0.0057 of FHP. 
The error arises from the luminosity density (±15%), from M/L (±20%) and 
from IMF (±20%), which are added in quadrature. 

3. Metal abundance in stars 

In the analysis of Kauffmann et al. (2003), the metallicity is also an output, but 
the results are not available. We here use the oxygen abundance determined from 
HII regions for nearby galaxies compiled by Kobulnicky & Zaritsky (1999). The 
metallicity shows a correlation with the luminosity of galaxies. Using oxygen as 
the metallicity indicator, integration over the LF in the B band yields, 

Z = lO 6 - 6±ai5 M Mpc" 3 , (2) 

where the solar composition is assumed. The zero point is set by the solar 
value, log[0/H]+12 = 8.83 at (Z/X) Q =0.0230 (Grevesse & Sauvel 2000), or 
Z & = 0.0163 using X = 0.71. 

There is an additional storage of heavy elements in white dwarfs, which 
are liberated only by Type la supernovae. From the abundance of white dwarfs 
(Bahcall & Soneira 1980), we estimate the heavy element (C+O) abundance as 

Z = 10 7 - 6 M o Mpc- 3 , (3) 



APS Conf. Ser. Style 



3 



which is much larger than (2). Heavy elements frozen in neutron stars are also 
a large amount s=s 1O 6 - 9 M Mpc~ 3 . 

4. Neutral and molecular gas mass 

FHP adopted the HI observation of optically selected galaxies by Rao and Briggs 
(1993), which yielded fi H i = (2.1±0.6) x 1CT 4 at h = 0.72. Since then, a number 
of blind HI surveys were carried out. Among them the largest sample (1000 
galaxies) was obtained by the HIPASS survey (Zwaan et al. 2003), which gives 
^hi = (4.2 ± 0.7) x 10~ 4 , twice higher than the value of Rao & Briggs. With 
the correction for helium (for both HI and H2), the amount of atomic gas is 

^m+Hci = (6.2 ±1.0) x 10" 4 . (4) 

The molecular hydrogen abundance is estimated from the CO survey of 
Keres, Yun & Young (2003): 

n U2 = 1.6 ±0.6 x 10 -4 . (5) 

This is compared to Qn 2 = 2.1 ± 0.6 x 10~ 4 (FHP) obtained by summing the 
mean H2 abundance for each morphological class of galaxies (Young & Scoville 
1991) weighted by the abundance of morphologically classified galaxies. 

5. Hot gas in clusters 

In FHP the hot gas abundance in clusters was estimated by integrating the 
cluster abundance for mass M > 1 x 10 14 /iM Q (Bahcall & Cen 1993) and mul- 
tiplying the gas fraction obtained from X ray observations. The cluster mass 
was defined by the Abell radius. Now, the advancement in cluster studies allows 
us to use the mass within r < r2oo, where r2oo is the radius at which matter 
density p = 200p cr it. From a theoretical ground this may give a better measure 
for the mass of the virialised system. Reiprich & Bohringer (2002) estimated 
from ROSAT All-Sky Survey that Q d = °- 012 ^ooo4 for clusters with mass larger 
than M = 4.5 x 1O 13 M0, which are visible with X rays. 

The cosmic value of the baryon to total mass ratio from the WMAP (Spergel 
et al. 2003) is 

n b /Q m = 0.178(1 ±0.09). (6) 

We estimate the ratio of the stellar to total mass from the mean value of Mj Lb = 
(450±100)/i and M S /L B = 4.5(1±0.20), giving M s /M tot = 0.014(1±0.30). This 
is somewhat larger than the stellar mass density (1) divided by the total matter 
density from WMAP Q rn = 0.26 ± 0.05: ft s tar/^ m = 0.010 ± 0.004. Assuming 
that the baryon to dark matter ratio in clusters agrees with the cosmic value 
and subtracting the stellar mass from the total baryonic mass, we estimate the 
hot gas abundance: 

fid gas = 0.0020 ± 0.0006. (7) 

The significant downward shift compared with FHP is due to the different defi- 
nition of the radius with which the cluster mass is defined. 
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6. Warm and cool plasma 

FHP inferred the presence of copious warm and cool plasma based on the univer- 
sality of the baryon to dark matter ratio at large scales, and suggested that this 
component fills the gap between the cosmic baryon abundance from Big Bang 
nucleosynthesis and that estimated from observed baryons in the local Universe. 
The evidence for abundant warm gas around galaxies was presented by the de- 
tection of O VI absorption in the UV spectrum (Tripp, Savage & Jenkins 2000). 
This year, WMAP gave an accurate estimate for the baryon abundance, which 
agrees with the value from the deuterium and helium abundance in Big Bang 
nucleosynthesis. This erases any doubts concerning the estimate of the cosmic 
baryon abundance from the nucleosynthesis argument. 

We may estimate the abundance of baryons associated with galaxies from 
the the mean M/L ratio and the LF of galaxies, assuming that the baryon to 
dark matter ratio is universal when averaged over large scales. The M/L of 
Milky Way is known to be ~100 at 200 kpc (e.g., Kuijken 2003). The analysis 
of Prada et al. (2003) (see also Zaritsky et al. 1997) using the motion of 3000 
satellites around host galaxies derived from the SDSS yielded (M/L) = 120/i 
at the 'virial radius'. The least model-dependent method to measure the mass 
associated with galaxies is to use gravitational lensing shear around those with 
known redshifts. McKay et al. (2001) estimated the galaxy mass using the SDSS 
sample, giving (M/L r ) = (170 ± 21)/i for R < 260 kpc from the r band data. 
(The i band data give a smaller value, and g band data give a larger value.) 

These M/L values are significantly (by about a factor of 2) smaller than 
those for clusters. Taking the lensing value of M/L r = (170 ± 20)h and the r 
band luminosity density, we estimate Vt m = 0.14 ±0.02 for the matter associated 
with galaxies (within the virial radius). This leads to Qb = 0.025 when multiplied 
by the universal value of (6). Subtraction of O s and f2Hi+Hei+H 2 gives 

^w/c gas = 0.022 ± 0.005 (8) 

for the warm baryon component around galaxies. 

An alternative path to estimate the warm/cool baryon abundance is to 
subtract stars, neutral and hot ionised gas from the global baryon amount, 
which is accurately known after the WMAP observation: 

^w/c gas = 0.044 - 0.0025 - 0.0020 - 0.0008 = 0.039 ± 0.004. (9) 

The discrepancy of (8) and (9) implies 'missing baryons': the gap between 
the two estimates suggests the presence of baryons that are not immediately 
associated with galaxies. We do not count in (9) cool (sa 10 4 K) baryons in 
Lyman a clouds, which were estimated to give ft = 0.002 ± 0.001 (FHP), but 
this contribution is much too small to fill the gap. These missing baryons may 
be in the vicinity of galaxies beyond a few hundreds of kpc, or associated with 
dark clumps which do not shine as galaxies, as they occur in CDM simulations 
(Ostriker et al. 2003). The two possibilities may not be necessarily exclusive to 
each other. A possibility is not excluded that the missing baryons are present 
as a highly ionised diffuse component, though this is not very likely (see below). 

There is also a gap in the dark matter abundance between the cosmic value 
fidm = 0.22 and the amount associated with galaxies fidm ~ 0.12. According to 
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the hierarchical clustering calculation, we anticipate 15% of baryons are unbound 
(M < 1O 6 M ) and 25% are in clumps of mass 1O 6 M < M < 1O 1O M . CDM 
simulations (Ostriker et al. 2003) predict a lot of dark, low mass clumps. 

It would be interesting to note that if we adopt M/L r ~ 320/i (M/Lb — 
450/i) of clusters as the universal M/L of galaxies to estimate the global mass 
density, we would obtain £l m ~ 0.27 which is consistent with the cosmic mass 
density from WMAP. This implies that dark clumps are integrated into clusters, 
leading to a large M/L ratio, compared to that for galaxies, for which those dark 
components are excluded from accounting. This suggests that the majority of 
dark clumps reside in outskirts of galaxies, or in filaments and groups of galaxies, 
so that dark matter component that is not counted in the estimate of C * (M/L) 
is localised, rather than smoothly distributed through the Universe, and so are 
baryons. 

7. Metal abundance outside stars 

Taking the metallicity of interstellar gas given in sect. 2, we find Z ~ 1.5 x 
1O 6 M . For the cluster gas, we infer Z m 1.6 x 1O 6 M , adopting 1/3 solar. 
Very little is known for warm and hot plasma. If we assume the heavy element 
abundance of 0.01 solar as in globular clusters, or in typical Lyman a clouds, 
we get Z « (0.5 — 0.9) x 1O 6 M . Therefore, the metal abundance is dominated 
by that in dead stars in galaxies. 



Table 1. Summary of the cosmic baryon budget 



component 


FHP 


new estimate 


stars 


0.0019-0.0057 


0.0025±0.0008 


HI+Hel gas 


0.00025-0.00041 


0.00062±0.00010 


H2 molecular gas 


0.00023-0.00037 


0.00016±0.00006 


hot plasma in clusters 


0.0014-0.0044 


0.0020±0.0006 


warm and cold plasma (by sum) 


0.0072-0.030 


0.022±0.005 


(by subtr.) 




0.037±0.004 


total 


0.011-0.041 


0.044±0.004 



8. Summary 

In this report I have presented an updated accounting of the cosmic baryon 
budget, and have given an estimate for the heavy element abundance. The 
summary of the cosmic baryon budget is shown in Table 1. Each entry does 
not differ greatly from that given in FHP, but an increased accuracy reveals 
a missing baryon component which amounts to (0.37 — 0.22)/0.44 « 35% of 
the total. This would provide an interesting test for cosmic simulations of the 
structure formation. 
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